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Abstract 

The Gribov copies and their consequences on the infrared behavior of the gluon 
propagator are investigated in Euclidean Yang- Mills theories quantized in linear co- 
variant gauges. Considering small values of the gauge parameter, it turns out that 
the transverse component of the gluon propagator is suppressed, while its longi- 
tudinal part is left unchanged. A Green function, Gtr(k), which displays infrared 
enhancement and which reduces to the ghost propagator in the Landau gauge is 
identified. The inclusion of the dimension two gluon condensate (AV) is also con- 
sidered. In this case, the transverse component of the gluon propagator and the 
Green function Gtr(k) remain suppressed and enhanced, respectively. Moreover, the 
longitudinal part of the gluon propagator becomes suppressed. A comparison with 
the results obtained from the studies of the Schwinger-Dyson equations and from 
lattice simulations is provided. 
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1 Introduction 



Gribov ambiguities [XJ and their relevance for the nonperturbative aspects of Euclidean 
Yang-Mills theories have witnessed growing interest in recent years. These ambiguities, 
affecting the Faddeev-Popov quantization formula, deeply modify the infrared behavior 
of Yang-Mills theories, as one learns from the large amount of results obtained in the 
Landau gauge El El El El E| as well as in the Coulomb gauge [HI El CHI El EE] • 

As pointed out in |T3|, the existence of these ambiguities is due to the lack of a globally well 
defined gauge-fixing procedure. Among the class of covariant gauges, the Gribov ambigui- 
ties have been much investigated in the Landau gauge, where the gauge field is transverse, 
d^A® = 0. This property plays an important role here. It ensures that the Faddeev-Popov 

operator, Ai ab (A) = —d^ ( K d^5 ab — gf abc A c ^j, is hermitian, M. = Ai^ . Its eigenvalues are 
thus real. Concerning now the quantization of Yang-Mills theories in the Landau gauge, 
it turns out that, as a consequence of the existence of Gribov copies, the domain of inte- 
gration in the Feynman path integral has to be restricted to the so called Gribov region Q 
[H El El HI El El Ej, which is the set of all transverse fields for which the Faddeev-Popov oper- 
ator is positive definite, i.e. Q = {A^, d^A* = 0, M ab {A) = -<9 M (d^ - gf abc A c ^ > o}. 
The boundary dfl, where the first vanishing eigenvalue of the Faddeev-Popov operator ap- 
pears, is called the first Gribov horizon. For the partition function of Yang-Mills theories 
in the Landau gauge, one has 

Z = [ DA5(dA) (det(-dH ab + gf abc Ald^)e^5 d * xFa ^ 

= / DADcDc5(M)e-^ d4 ^ F - + / d4 ^ a ^( 9 ^ a6 - 9/Qbc ^) cb . (1.1) 
Jq 

The restriction of the domain of integration to the region Q has important consequences on 
the infrared behavior of the gluon and ghost propagators |H El El HI El Ej • More precisely, 
in the Landau gauge, the gluon propagator (^A^tyA^— k)^ is found to be suppressed in 

the infrared, while the ghost propagator (c a (k)c b (—k)) is enhanced, according to 



(Al(k)Al(-k)) = *- (V - 4£) ^ , (1.2) 



and 



<W*) = J^J2^ b (c a (k)c b (-k) 



Qgh(k) k ^0 ~ ^4 • (1-3) 

The Gribov parameter 7 in eq. ([1.2)1 has the dimension of a mass and is determined by a 
gap equation which, to the first order approximation, reads 

3Ng2 f d ' k 1 ! . (L4 ) 



4 J (2vr) 4 fc 4 + 7 4 
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From equation (jl.3|) . one sees that the infrared behavior of the ghost propagator is more 
singular than the perturbative behavior 1/k 2 . This infrared enhancement is known as the 
Gribov-Zwanziger horizon condition [HIUE], generally stated as linu^o {k 2 Q g h{k)) = 0. 
Remarkably, this behavior of the gluon and ghost propagators in the Landau gauge has 
received many confirmations from lattice simulations O E3 HU H3 [TBI EH I2H I2H 122] 
as well as from recent studies of the Schwinger-Dyson equations [2l?l 12^ 121^ I27| I2l?l I2T?1 13TJ] . 

Several results have been established about the Gribov region Q. It has been proven 
that Q is convex and bounded in every direction [21 E|. Moreover, every gauge orbit 
passes inside Q [E]. The latter result is deeply related to the possibility of introducing an 
auxiliary functional, J 7 {A) = J ^xA^A^, whose minimization with respect to the gauge 
transformations provides a characterization of the region fl [2JEII- It can be checked in 
fact that the Gribov region fl can be identified as the set of all relative minima of the 
functional J 1 {A) . It should also be mentioned here that the region fl is not free from Gri- 
bov copies, i.e. Gribov copies still exist inside fl (HI EH E2] • To avoid the presence of these 
additional copies, a further restriction to a smaller region A, known as the fundamental 
modular region, should be implemented jBJ EH1 E2] • The region A is contained in the Gri- 
bov region fl, being defined as the set of all absolute minima of the auxiliary functional 
F{A). However, it is difficult to give an explicit description of A. Recently, it has been 
argued in [33] that the additional copies existing inside fl have no influence on the expec- 
tation values, so that averages calculated over A or fl are expected to give the same result. 

Besides the Landau gauge, the effects of the Gribov copies on the infrared behavior 
of Yang-Mills theories have been studied to a great extent in the noncovariant Coulomb 
gauge [3 El CHI HI1 U2] , diA°; — 0,i = 1,2,3. In particular, the Coulomb gauge allows for 
a direct study of the confining properties of the static potential V(r) between an external 
quark pair at spatial separation r. It turns out that V(r) can be accessed by means of 
the 44-component, (^A%(l&,t)A\(0)\ of the gluon propagator. Moreover, in analogy with 

the Landau gauge, the spatial components of the gluon propagator, (^A1{k)A h -{—k)^, are 
found to be suppressed in the infrared, a feature confirmed by lattice simulations |3H EE] • 

Concerning now other covariant gauges, although the presence of the Gribov copies is 
certainly to be expected [T3], as explicitly confirmed by a recent study of the zero modes 
of the Faddeev-Popov operator in the maximal Abelian gauge very little is known 
about their consequences on the infrared behavior of the gluon and ghost propagators 
and, more generally, on the whole Yang-Mills theory. The need for such an investigation 
is motivated by the increasing belief that the Gribov copies might have a crucial role for 
the infrared region of Yang-Mills theories as well as for color confinement. It would be de- 
sirable thus to improve our understanding on how the Gribov copies manifest themselves 
in different gauges and how they modify the infrared behavior of the theory. 

This work aims at studying this issue in the class of covariant linear gauges, d^A^ = —ab a , 
where a is the gauge parameter and b a is the Lagrange multiplier. The task is far from 
being trivial. Many features of the Landau gauge are lost for a generic nonvanishing a. 
The Faddeev-Popov operator, Ai ab (A) = — <9 M [p^5 ab — gf abc A c ^, is now not hermitian. 
Moreover, a suitable minimizing functional in these gauges is not yet at our disposal. 
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As a consequence, the identification of the region to which the domain of integration in 
the path-integral should be restricted becomes difficult to be handled. All this makes a 
complete treatment of the Gribov copies in linear covariant gauges far beyond our present 
capabilities. Nevertheless, we shall be able to establish some preliminary results which 
enable us to obtain a characterization of the infrared behavior of the gluon propagator, 
at least for small values of the gauge parameter a. Considering in fact small values of the 
gauge parameter a, will allow us to stay close to the Landau gauge fixing. The present 
covariant gauge can be considered thus as a kind of deformation of the Landau gauge, 
whose results have to be recovered in the limit a — > 0, thereby providing a useful consis- 
tency check. 

The output of our findings can be summarized as follows. As in the case of the Lan- 
dau gauge, the transverse component of the gluon propagator turns out to be suppressed 
in the infrared. Moreover, its longitudinal part remains unchanged. 

Concerning now the behavior of the ghost fields in linear covariant gauges, it turns out 
that, instead of the ghost propagator, the Green function which is enhanced in the infrared 
is given by the quantity Gtr(k), defined as 



Gtr(k) 



1 



N 2 



where 



.4 



aT 



M-\A 7 



<9/A 
d 2 



ab 



Ad 
A V 1 



;i.5) 



1.6) 



is the transverse component of the gauge field and LA4 1 ( A T ] 
the Faddeev- Popov operator M. ab (A T ), 



ab 



M ab (A T ) = -dp (<V ab - 9f abcA 



stands for the inverse of 



;i.t) 



The Green function Qtr{k) is found to obey the Gribov-Zwanziger horizon condition, i.e. 
linifc^o (k 2 Qtr(k)) 1 = 0. It should be remarked that Gtr(k) does not coincide with the 
ghost propagator for a generic value of the gauge parameter a. However, Qtr(k) reduces 
to the ghost two-point function for vanishing a, so that our results turn out to coincide 
with those of the Landau gauge in the limit a — ► 0. 



The infrared behavior of the gluon propagator and of Qtr{k) will be investigated also 
in the presence of the dimension two gluon condensate (ApAfy. A detailed study of this 

condensate in linear covariant gauges can be found in [37]. In the presence of (A^A^j, 
the infrared suppression of the transverse component of the gluon propagator is enforced. 
Furthermore, its longitudinal component turns out to be suppressed as well. Concerning 
the Green function Q tr (k), its infrared enhancement is not modified by the condensate 

Aa Aa K 



The work is organized as follows. In Sect. 2 a few properties of the Gribov copies in 
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the linear covariant gauges, and for small values of the gauge parameter a, are presented. 
In Sect. 3 the infrared behavior of the gluon propagator, of the Green function Qtr{k) and 
of the ghost propagator is worked out. Sect. 4 accounts for the inclusion of the gluon 
condensate (A a ll A a \ Sect. 5 is devoted to a comparison of our results with those obtained 
from the analysis of the Schwinger- Dyson equations and from lattice simulations. In 
Sect. 6 we present our conclusion. 



2 A few properties of the Gribov copies in the linear 
covariant gauges 

In this section we shall discuss a few properties of the Gribov copies in the linear covariant 
gauges. Let us begin by considering the expression of the Euclidean Yang-Mills action 
quantized in the linear covariant a— gauges, namely 

Sym + S 9f = \j d 4 xF* u F« u - J d 4 x (b a 8A a + ^b a b a + c a d, (d,5 ab - gf abc A 

(2.8) 

where c a , c a stand for the Faddeev-Popov ghosts. The Lagrange multiplier b a has been 
introduced to implement the gauge condition 

d„Al = -ab a , (2.9) 

from which it follows 

J d A x (b a dA a + ^b a b a ^j =>- ~ Jd A x (<V^) 2 • (2- 10 ) 

From the relation ()2.9|) we see that the field A a is not transverse, due to the presence of 
the gauge parameter a. Of course, in the limit a — > 0, we recover the Landau gauge, 
d^A® = 0. In what follows, it will be useful to decompose the gauge field A^ into transverse 
and longitudinal components, according to 



A$ = Af + Af, (2.11) 



with 



A aT ( r _ <VM A a 

H - Q2 J A U 7 

d,Af = 0, (2.12) 

and 



^ = dj ^K = -*%b« 



d»Af = -ab a . (2.13) 

As already underlined, we shall restrict ourselves to the case in which a is small, i.e. 
a <C 1, so that the longitudinal component A^ L in eq. f|2.11|) can be considered as a small 
perturbation with respect to the transverse part A a J . 



Let us proceed by proving the following statement. 
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Statement 

If the transverse component A a ^ of the gauge field = (yA a ^ + A^f^j belongs to the 

Gribov region ft, i.e. Af G ft, ft = [Af , d„Af = 0, - <9 M (<9 M <5 afe - gf abc Af) > o}, 

then the Faddeev- Popov operator Ai ab (A) = —d^ (d lx S ab — gf abc A c ^j has no zero 
modes. 

Proof 

The proof of this statement is done by assuming the converse. Suppose indeed that 
the operator Ai ab (A) has a zero mode, i.e. it exists a (p a (x,a) such that 

- d, (d,5 ab - gf abc Af - gf abc Af) ^ b (x, a) = , (2.14) 

which, from eq. f!2.13|) . becomes 

- d, [d,5 ab - gf abc Af + agf abc {^b^ ) $ b (x, a) = . (2.15) 

Since a is small we perform a Taylor expansion of (p a (x, a), namely 

oo 

i?(x,a) = y E l a n ift{x), (2.16) 

n=0 

with 

(p a o(x) = a (x,a)\ a=o . (2.17) 
Equation (|2.15jl splits thus in a tower of equations, given by 



- d, (d,5 ab - gf abc Af) (pl{x) = , (2.18) 



0„ 



- d, (d,5 ab - gf abc Af) (p b (x) - d, gf abc -f 6 C ifi(x) = , (2.19) 



d, (d^ - gf abc Af) (p b (x) - d, gf abc ( ^b c ) (p\(x) ) = , (2.20) 



d, (d,8 ab - gf abc Af) (p b (x) - d, gf abc ( -£b c ) (p b (x) ) = , (2.21) 



Q2 
Q2- 



and so on. 



Moreover, since A a ^ belongs to the Gribov region ft, A°^ G ft, it follows that ^Pq{x) 
in the first equation ()2.18|) necessarily vanishes, i.e. ^Pq(x) = 0, since the operator 
— (dfj,5 ab — gf abc A c ^ has no zero modes in ft. Furthermore, setting <Pq(x) = in 
the second equation (|2.19|) . we get 

- d, (d,5 ab - gf abc Af) (p\{x) = , (2.22) 

which, in turn, implies the vanishing of ^\{x), i.e. (f b (x) = 0. As a consequence, 
eq.flTHJ) reads 

- d, (d,5 ab - gf abc Af) (p\{x) = , (2.23) 



from which (p\{x) = 0. It is apparent thus that the argument can be easily iterated 
to the whole tower of equations, implying that (p b (x, a) vanishes, <p b (x, a) = 0. This 
concludes the proof of the statement. 
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In particular, from this result it follows that if belongs to the Gribov region, A a ^ G Q, 

the field = [Af + A° L ) does not possess Gribov copies which are closely related, i.e. 
obtained from through an infinitesimal gauge transformation 

Al = Al + (d,5 ab - gf abc Al) u b , (2.24) 

where u a (x) denotes the parameter of the gauge transformation. Indeed, from the condi- 
tion 

d,Al = d»A« , (2.25) 

we should have 

d, (d,5 ab - gf abc A^) u b = , (2.26) 
which, however, has no solution for uj a {x). 



2.1 Restriction of the domain of integration in the path- integral 

The results obtained in the previous sections suggest to restrict the domain of integration 
in the path integral to the region Q defined as 

fi = {A% Al = A? + Af, A? G n} , (2.27) 

i. e. Q is the set of all connections whose transverse part A a ^ belongs to the Gribov region 

Q = [Af, d^Af = 0, -dp (d^5 ab - gf ahc Af) > o}. This choice is motivated by the 
following considerations: 

• Since the gauge parameter a is small, a < 1, the region Q can be regarded as a 
deformation of the Gribov region Q. It is apparent in fact that in the limit a — > 
the region Q is recovered, namely 

lim O = n . (2.28) 

• As discussed before, the Faddeev-Popov operator, M ab (A) = — <9 M (d fl 5 ab — gf abc A c ^j , 

has no zero modes within Q. Therefore, the restriction to Q allows us to get rid of the 
Gribov copies which can be obtained through infinitesimal gauge transformations. 

• The effective implementation in the path-integral of the restriction of the domain 
of integration to the region Q can be done by repeating the no-pole prescription of 
Gribov's original work pQ . Indeed, since the transverse component A ^ of any field 
belonging to Q lies in the Gribov region Q, to implement the restriction to Q it 
will be sufficient to require that the Green function Qt r {k) of eq. (jl.5J) has no poles 
for any A a ^ G f2, except for a singularity at k 2 = 0, corresponding to the Gribov 
horizon dQ. 
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Thus, for the partition function of Yang-Mills theories in linear covariant gauges, we write 

Z = J IMDW(0A + aft) det ^ 

= J DA det (-d, (8 a % -gf abc A^)) e-tif* x V^ + &f* x ( dAa ?)v(n) , (2.29) 

where the factor V(Cl) implements the restriction to the region Q. An explicit character- 
ization of V(0) and of its consequences on the infrared behavior of the gluon propagator 
and of Gtr(k) will be discussed in the next section. 

Finally, it is worth to spend a few words on the aspects which remain uncovered by 
the present investigation. Even if the restriction to the region Q allows us to get rid of the 
Gribov copies which are closely related, i.e. attainable by infinitesimal gauge transfor- 
mations, we still lack a treatment of the copies which cannot be attained by infinitesimal 
transformations. This task is beyond our present possibilities, as the knowledge of a 
suitable auxiliary functional associated to the linear covariant gauges would be required. 
Nevertheless, since we are limiting ourselves to small values of a, the restriction to the 
region Q looks quite natural. 



2.2 Characterization of the factor V(f2) 

As already remarked, the factor V(Cl) can be accommodated for by requiring that the 
Green function Qtr{k) of eq. ()1.5|) has no poles for a given nonvanishing value of the mo- 
mentum k. This condition can be understood by recalling that the region Q is defined 
as the set of all transverse gauge connections |A^ a |, <9 M v4^ a = 0, for which the operator 

M ab (A T ) is positive definite, i.e. M ab (A T ) = -<9 M (d^ ah - gf abc A T ^ > 0. This implies 

that the inverse operator —d^ (d^5 ab — gf abc A^ c ^j , and thus Qtr(k), can become large 
only when approaching the horizon dQ, which corresponds in fact to k = p. The 
quantity Qt r (k) can be evaluated order by order in perturbation theory. Repeating the 
same calculation of [1 j, we find that, up to the second order 



fc <*> k swiU • (2 - 3o) 



with p(k, A T ) given by 



N 9 2 kfik u 1 A Ta/ \ A Ta 



m _ lv „ (2.31) 

V being the Euclidean space-time volume. According to PP, the no-pole condition for 
Qtr(k) reads 



p(0,A T ) < 1 
1 

4N 2 -IV ^ q 2 



p^ aT ) = 7^CE4Ul a (gMr(-?)) • (2.32) 
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Therefore, for the factor V(Q) in eq. (|2.29|) we have 

V(n) = 6(l-p(0,A T )), (2.33) 
where 8(x) stands for the step function*. 

3 The gluon propagator 

In order to study the gluon propagator, it is sufficient to retain only the quadratic terms in 
expression (J2.29|) which contribute to the two-point correlation function (^A a ^(k) A h v (— k) 
Making use of the integral representation for the step function 



6(1 - p(0, A )) = / _U_ e v(i-P^ )) (3>34) 
J-ioo+e 2mr) 

for the partition function ()2.29|) one get 

£ qU adr = N J D A^-e^ Xo ^ e< s ^ +vp{a ' AT) ) , (3.35) 

where Af is a constant factor and .Squadr stands for the quadratic part of the quantized 
Yang-Mills action 

S quadr = \ £ (q 2 A*(q)A^-q) - (l - i) ^g^(g)^(-g)) . (3.36) 

From 

Al a {q)Al a {-q) = ^ - ^)(<W ~ q ~f] Al{q)Al{- q ) 

= A«(q)A«(-q) - ^A^AH-q) , (3.37) 

it follows that 



iW = AT I Z)^e^ b "e-^/J (9)2 -^ ( - 9) , (3.38) 
J 27fi 



drj 

^quadr = J\ J 

with 

* - {h i^M) «- - - fr - 3 + «- • <"»> 

Integrating over the gauge field, one has 

iw =x] ^~ iogv ( det Qty h =^'j ^ m . ( 3 - 4 °) 

where /(r/) is given by 

/(^) = r/-logr / -^(AT 2 -l)Y:iog^ 4 + ] ^^ . (3.41) 



*0(a;) = 1 if x > 0, 0(x) = if x < 0. 
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Following P, the expression (|3.40|) can be now evaluated at the saddle point, namely 

^quadr « , (3.42) 

where r] is determined by the minimum condition 

1 3 No 2 „ 1 

/0 9 + TV 2 — 1 2V J 

Taking the thermodynamic limit, V — > oo, and introducing the Gribov parameter 7 pQ 

4 VoNg 2 1 , , 

7 = — n , K — > 00 , 3.44 

- 1 2V 



- N 9 2 I tAi-T-TZ = 1 > ( 3 - 45 ) 



we get the gap equation 

3 2 /• A 1 
4^7 (2vr)V + 7 4 

where the term 1/770 in (|3.43|) has been neglected in the thermodynamic limit. To obtain 
the gauge propagator, we can now go back to the expression for Z qua dr which, after 
substituting the saddle point value 77 = 770, becomes 

Z quadl = M J DAe^ AimtuAlt-,) ^ {3M) 



with 



Evaluating the inverse of Q'z, , we get the gluon propagator 

A°M)Al(- q) ) = ^ ( V - S^) + ^) . (3.48) 

A few remarks are now in order. 



Let us first observe that the gap equation ()3.45|) defining the parameter 7 has the same 
form of that obtained by Gribov in the Landau gauge pQ. This is an expected result, since 
the factor p(0, A T ) in equation (J2.32j) depends only on the transverse component A^ a . 

As it is apparent from the expression ()3.48|) . the transverse component of the gluon prop- 
agator is suppressed in the infrared, while the longitudinal component is left unchanged. 
Moreover, as we shall see later, the behavior of the longitudinal part turns out to be 
modified once the gluon condensate (A^A^ is taken into account. 

Finally, in the limit a — > 0, the results of the Landau gauge are recovered. 
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3.1 The infrared behavior of Qt r {k) 

Let us discuss now the infrared behavior of the Green function Qtr{k) of eq. (jl.5j) . which 
is obtained upon contraction of the gauge fields in the expression ()2.3H) . namely 

fc( *> 53 hr^m ■ (3 - 49) 



with 



P(k) = 9 2 ^-r^ i ,f* 4 „ 1 , 2 (Al\q)Al\- q )) . (3.50) 



iV 2 -l k 2 J (2n) 4 (k-q) 
From the expression of the gluon propagator in eq. (j3.48J) , it follows 

Mpk v f d 4 q 1 q 2 ( q^q v 



P(k ) = g N— J — T U - . (3.51) 



(2jr)*(fc- g )V + 7 4 V * 9 2 
Furthermore, from the gap equation (|3.45|h it turns out 



d 4 q 1 ( s q^q 
(2tt)V + 7 



N 9 2 I T^A-l- -i <W - ^ = <V > ( 3 - 52 ) 



so that 

>k^k v f d 4 q k 2 -2qk 1 / q^ 



i-M'^J^W^Tt*^- 1 ?- (3 ' 53) 



Note that the integral in the right hand side of eq. fl3.53j) is convergent and nonsingular at 
k = 0. Therefore, for k ~ 0, 

(1-#)U-^ 2 . (3-54) 
where X stands for the value of the integral 

d 4 q 1 1 1 

(2vr)V(<? 4 + 7 4 ) ~ 327f 

Finally, for the infrared behavior of the Green function Qt r (k) we get 



Z ~ I /« \4„2(Mj.nA\ - <?9^2 • ( 3 - 55 ) 



«M« « 5a^iF • (3-56) 

One sees thus that Qtr{k) is enhanced in the infrared, obeying the Gribov-Zwanziger 
condition lim fc ^ {k 2 Qtr{k)) 1 = 0. 



3.2 Analysis of the ghost propagator 

For the sake of completeness, let us discuss here the infrared behavior of the ghost two- 
point function, G g h{k), given by 
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with 

Making use of the expression for the gluon propagator in eq. ([3.48)1 and of the equation 
(13.52)1 . it follows that, in the infrared, 

3Ng 2 k 2 aNg 2 f d 4 q kjjk - q) v q^q v 
1 " W 128^ " ^ 7 (^T^T^T • (3-59) 



The second term in the right hand-side of eq. fl3.59j) can be easily evaluated by means of 
the dimensional regularization. Adopting the MS scheme, the final expression for the 
factor (1 — u}(k)) is found 



3Ng 2 k 2 aNg 2 , k 2 

uik)) « — V log^ . 3.60 

128vr 7 2 64vr 2 y /Z 2 



One sees that, in the present case, the Gribov-Zwanziger horizon condition is jeopardized 
by the term log=z in expression (|3.6U|) . As it is apparent from the presence of the gauge 
parameter a, this term is due to the contribution of the longitudinal components of the 
gauge field. Note that the longitudinal modes do not contribute to the Green function 
g tr (k) in eq.flUHD. 



Inclusion of the dimension two condensate (A^A^ 



The dimension two gluon condensate {A^A^j has received much attention in the last 

years |SHl EHl E01 SH E21 SH SSI SHI SZ1 SHI- This condensate turns out to contribute 
to the gluon two-point function, as observed in jl^l within the operator product expan- 
sion. As such, it has to be taken into account when discussing the gluon propagator. A 
renormalizable effective potential for I^A a ^A a ^ in linear covariant gauges has been con- 
structed and evaluated in analytic form in [37|. The output of this investigation is that 
a nonvanishing value of the condensate (A^A^J is favoured since it lowers the vacuum 
energy. As a consequence, a dynamical tree level gluon mass is generated in the gauge 
fixed Lagrangian [H7j. The inclusion of the condensate (^A^Af^) in the present framework 
can be done along the lines outlined in [20], where the effects of the Gribov copies on 
the gluon and ghost propagators in the presence of (A a ^A a ^) have been worked out in the 
Landau gauge. Let us begin by giving a brief account of the dynamical mass generation 
in linear covariant gauges. Following jSTj, the dynamical mass generation in these gauges 
is described by the following action 

S(A,a) = S YM + S gf + S <r , (4.61) 

where Sym, S g f are the Yang-Mills and the gauge fixing terms, as given in eq. f|2.8|) . The 
term S a in eq. (14.61)1 contains the auxiliary scalar field a and reads 
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The introduction of the auxiliary field a allows us to study the condensation of the local 
operator A^A^. In fact, as shown in [37j, the following relation holds 

{°) = -\{KK) ■ ( 4 - 63 ) 

The dimensionless parameter ( in expression (J4.62)) is needed to account for the ultraviolet 
divergences present in the vacuum correlation function (A 2 (x)A 2 (y)). For the details of the 
renormalizability properties of the local operator A^A^ in linear covariant gauges we refer 
to jSTj. The action S(A, a) is the starting point for evaluating the renormalizable effective 
potential V(a) for the auxiliary field a, obeying the renormalization group equations. The 
minimum of V(cr) occurs for a nonvanishing vacuum expectation value of the auxiliary 
field [HZ], 2.e. (a) ^ 0. In particular, from expression (|4.61j) . the first order induced 
dynamical gluon mass is found to be 

m 2 = , (4-64) 

Co 

where Co is the first contribution to the parameter (, given by [37] 

C = ^ + Ci + o(^ 2 ), 
g 2 



3 (78 - 26a 2 + 3a 3 + 18a log a) (N 2 - 1) 
C ° " 2(3a-13) 2 N^' (U) ' J) 



We remind here that in linear covariant gauges, the Faddeev-Popov ghosts c a , c a remain 
massless, due to the absence of mixing between the composite operators A a ^(x)A a ^(x) and 
c a {x)c a (x). This stems from additional Ward identities present in these gauges [5T] . which 
forbid the appearance of the term c a 



4.1 Infrared behavior of the gluon propagator in the presence 

of {A*A$ 

It is worth underlining that the action S(A, a) leads to a partition function which is 
still plagued by the Gribov copies. It might be useful to note in fact that the action 
(Sym + S a ) is left invariant by the local gauge transformations 

8AI = -Z^V, (4.66) 
5a = gA*d,oo\ 

6(S YM + S a ) = 0. (4.67) 

Therefore, implementing the restriction to the region Q, for the partition function we 
obtain 

Z = J DA Dadet (-8, (5 a % + gf ahc A$) e -(i/ d4 ^ i? -+^/^(^) 2 +s CT ) v( ^ ) ? 

(4.68) 

with the factor V(f2) given in eqs. ()2.32|) . (j2.33|) . To discuss the gluon propagator we pro- 
ceed as before and retain only the quadratic terms in expression ()4.68j) which contribute to 
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the two-point correlation function (^A^{k)A° v (— k)). Expanding around the nonvanishing 
vacuum expectation value of the auxiliary field, (a) ^ 0, one easily get 

= ( DA *L^*n e -\Y.,*%<M)<&*™ (4 . 69) 
J 2ni 

where the factor Q°^ v is now given by 



(4.70) 

-2W = AT / l^- 1 ^ (det Q J) "* =Af'J , (4.71) 



Integrating over the gauge field, one has 

^quadr J v / 

with 



/(^)=r ? -logr ? -^(iV 2 -l)^log^ 4 + mV + ] ^^^ 7 ) ■ (4.72) 
Evaluating -Z qua dr at the saddle point, yields 

^quadr « ^ , (4.73) 

where 770 is determined by the minimum condition 

i_I_^V 1 = (4 74) 

Taking the thermodynamic limit, V — > 00, and introducing the Gribov parameter 

we get the gap equation in the presence of the dynamical gluon mass, corresponding to a 
nonvanishing condensate (A^A 1 ^), namely 

d 4 q 1 



= 1 (4 76) 

(2vr) 4 g 4 + m 2 ? 2 + 7 4 ' v ' ; 

Note that the dynamical mass m appears now explicitly in the gap equation (J4.76|) . To 
obtain the gauge propagator, one goes back to the expression (|4.69|) which, when evaluated 
at the saddle point value r\ = i] , yields 

iW =M j DAe-^ A W<&*™ , (4.77) 

with 

fl- = + m > + J) ( V - &*) + (£ + „,.)),- . (4.78) 

Thus, for the gauge propagator in the presence of the dynamical gluon mass m we get 

(A^AK-q)) = 5" b ( f 4 - M^) + • (4.79) 

x ^ ' + m^g^ + 7 4 y g^ / g^ + am 2 q 2 J 

We note that, due to the presence of the mass m, the infrared suppression of the transverse 
component of the gluon propagator is enforced. Moreover, also the longitudinal component 
gets suppressed. 
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4.2 The infrared behavior of Qtr{k) in the presence of (A^Az) 

It remains now to discuss the infrared behavior of the Green function Qtr{k) in the presence 
of (A^A^S. This can be easily worked out by repeating the analysis done in the previous 
sections. From the expression of the gluon propagator ([4.79)1 . it follows that 



with 



N 2 -l k 2 J (2n) (k - qY 
kykv f d 4 q 1 q 2 ( q^ 



k 2 J (2vr) 4 (k - qf q A + m 2 q 2 + ^ ^q 2 1 " ' LS1 ' 



Also, from the gap equation (|4.7fj)l . one has 

d A q 1 L _ QtJl, 

(2tt) 4 q 4 + m 2 q 2 + 7 4 V " 9 



N 9 2 I THI ^J,^ ( V - ^ ) = <W > (4-82) 



so that 

74 „ J.2 



i_ m - v /" g k 2 -2qk l / 

W - ^ p y (27r) 4 (A ,_ g) 2 g 4 +m 2 ? 2 +T 4^ q 2 



3Ng 2 J 
'k^o ~ 4 

where J7 stands for the value of the integral 



(4.83) 



Thus, for k ~ 0, 



7 - [ d4q 1 

i (27r)V(g 4 + mV + 7 4 ) ' (4 ' 85) 

which is ultraviolet finite. Therefore, for the Green function Q tr (k), we get 

«W*»~"sivjvF- (486) 

exhibiting the infrared enhancement which, thanks to the gap equation ()4.76|) . turns out 
to hold also in the presence of the gluon condensate ( A^A"^ 



5 Comparison with the results obtained from lattice 
simulations and from the Schwinger-Dyson equa- 
tions 

Having investigated the infrared behavior of the gluon propagator and of the Green func- 
tion (Q tr (k)), as summarized by equations ()4.79|) and (|4.86|) . it is useful to make a com- 
parison with the results already available from lattice simulations and from the studies of 
the Schwinger-Dysons equations. Let us begin with the lattice data 
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5.1 Comparison with the lattice data 



In a series of papers |52] 153] 151] . Giusti et al. have managed to put the linear covariant 
gauges on the lattice. This has allowed for a numerical investigation of the transverse 
as well as of the longitudinal component of the gluon propagator. Following [53] , let us 
introduce the transverse and longitudinal form factors -Dr(<?) and Di{q) through 

(A°MAH- q) ) = S* (^M _ S*) + £M^ . (5 .87) 

The results obtained in [52] |53] |53| show that both D T (q) and D L (q) are suppressed in 
the low momentum region, see for instance Fig.3 and Fig.4 of |54j . Our results are in 
qualitative agreement with the lattice data. Indeed, from the expression (]4.79|) . we obtain 

D T (q) 



c/ 4 + m A q J ' + 7 4 

DM = 2 °f a , (5-88) 
q l + am 2 

exhibiting infrared suppression. Note that, at least within the approximation considered 
in the present work, the suppression of the longitudinal form factor -Dl(c/) in eq. fl5.88j) 
is a consequence of the dynamical gluon mass, due to the gluon condensate (A®Afy, as 
already pointed out in [37] . Concerning now the ghost propagator and the Green function 
Qtr{k), to our knowledge, no results from lattice data are available so far. 



5.2 Comparison with the results obtained from the Schwinger- 
Dysons equations 

The infrared behavior of the gluon and ghost propagator has been investigated within 
the Schwinger-Dyson framework in [55] . Here, a power-law Ansatz for the transverse 
and longitudinal form factors of the gluon propagator as well as for the ghost form factor 
D gh (q) has been employed, according to 



D T (q) « (q 2 ' 
D L (q) « (q 2 Y , (5.89) 



and 

[c a (q)c b (-q) 



D gh {q) 



q 2 



D gh (q) « [q ) , (5-90) 

The results obtained for the infrared exponents (cr, p, (3) turn out to be similar to those 
of the Landau gauge, namely''' 

a > , 

P > o, 

= f = f • (5 - 91) 
^The explicit values of these infrared exponents as well as their dependence from the gauge parameter 
can be found in |55) . 
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indicating an infrared suppression of the transverse and longitudinal gluon form factors, 
and an infrared enhancement of the ghost propagator. Concerning the gluon propagator, 
these results are in qualitative agreement with our results as well as with the lattice 
data. However, concerning the ghost propagator, we have found that, instead of the 
ghost form factor, the quantity which is enhanced in the infrared is Q tr (k). For a better 
understanding of this point, it is worth reminding here that the result for the infrared 
exponents in eq. fl5.91|) has been obtained by using a bare- vertex truncation scheme Hoi. 
This approximation has been proven successful in the Landau gauge [23 123 123 123 EE 
I2H]- In particular, in the Landau gauge, no qualitative difference has been found if bare 
vertices are replaced by vertices dressed according to the Slanov- Taylor identities. This 
feature of the Landau gauge is believed to be deeply related to the nonrenormalization 
theorem of the ghost-antighost-gluon vertex, which holds to all orders of perturbation 
theory [5211211. Recently, the nonrenormalization theorem of the ghost-antighost-gluon 
vertex in the Landau gauge has been investigated through lattice simulations in 
which have provided indications of its validity beyond perturbation theory. However, to 
our knowledge, no such a theorem is available in linear covariant gauges, for a nonvanishing 
value of the gauge parameter a. Furthermore, according to the authors 55 , it is yet an 
open question whether the values of the infrared exponents in eq. (|5.9H1 remain unchanged 
if bare vertices are replaced by dressed ones. Our results suggest that a different behavior 
might be expected when dressed vertices would be employed. 



6 Conclusion 

In this work we have attempted at analyzing the effects of the Gribov copies on the gluon 
propagator in linear covariant gauges. By considering small values of the gauge parameter 
a, a few properties of the Gribov copies have been established, allowing us to investigate 
the infrared behavior of the gluon two-point function. 

As in the case of the Landau gauge, it turns out that the transverse component of the 
gluon propagator is suppressed in the infrared. Moreover, the longitudinal part is left un- 
changed, as shown in eq. (J3.48|) . The infrared behavior of the gluon propagator has been 
investigated also in the presence of the gluon condensate (^A a ^A a f S. In this case, the in- 
frared suppression of the transverse component is enforced. Furthermore, its longitudinal 
component turns out to be suppressed as well, as expressed by eq. (|4.79|) . These results 
are in qualitative agreement with those obtained from lattice simulations and from the 
analysis of the Schwinger-Dyson equations. 

Concerning now the behavior of the ghost fields in linear covariant gauges, the output of 
our analysis is that, instead of the ghost propagator, the Green function which exhibits 
infrared enhancement is given by Qt r (k), as defined in eq. ()1.5jl . It should be remarked 
that Qtr(k) does not coincide with the ghost propagator for a generic value of the gauge 
parameter a. However, Gtr{k) reduces to the ghost two-point function for vanishing a, 
so that our results turn out to coincide with those of the Landau gauge in the limit a — > 0. 

Needless to say, many aspects of the covariant linear gauges remain still to be inves- 
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tigated. A partial list of them is: 



As already pointed out, a suitable auxiliary functional corresponding to the linear co- 
variant gauge fixing condition is not yet at our disposal. As in the case of the Landau 
gauge El EH] , this functional could be very helpful for a characterization of the proper- 
ties of the Gribov copies not attainable by infinitesimal gauge transformations. 

From the present analysis, it emerges that the Green function Qtr{k) has a special role, 
as it obeys the Gribov-Zwanziger horizon condition and reduces to the ghost two-point 
function in the Landau gauge. Although its dependence from the transverse component 
A a ^ of the gauge field suggests that it might have a deeper meaning, it would be worth 
to have a better understanding of Qtrik). 

It would be useful to have a consistent framework to compute quantum corrections to 
the gluon propagator and to Q tr (k). This would amount to construct a local renormal- 
izable action in linear covariant gauges incorporating the effects of the Gribov copies, as 
done by Zwanziger in the Landau gauge jU [7j . 

Finally, it would be interesting to have more data in the linear covariant gauges from 
lattice simulations on the gluon and ghost propagators as well as on the Green function 
Qtrik). 
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